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SUMMARY 


This paper describes a new three-dimensional structured multiple-block volume grid generator called 
3DGRAPE/AL. It is a significantly improved version of the previously-released and widely-distributed program 
3DGRAPE, with many of the improvements taken from the grid-generator program 3DMAGGS 1 . It generates 
volume grids by iteratively solving the Poisson Equations in three-dimensions. The right-hand-side terms are 
designed so that user-specified grid cell heights and user-specified grid cell skewness near boundary surfaces 
result automatically, with little user intervention. Versatility was a high priority in this code's development, and 
as a result it can generate grids in almost any three-dimensional physical domain. Improvements include added 
kinds of forcing functions, improved control of cell skewness, improved initial conditions, convergence 
acceleration, the ability to take as input the output from GRIDGEN, and a simple but powerful graphical user 
interface(GUI). 


INTRODUCTION 


The original program, 3DGRAPE 2 *3, of which 3DGRAPE/AL is an updated version, is a batch-type 
program. It reads in pre-defined input data, generates the grid, and writes it out. For those boundary surfaces 
which are of interest (i.e., the body) it expects to read X,Y,Z coordinates of surface grid points which the user has 
pre-defined using other software. Other boundary surfaces of less interest (i.e., the outer boundary) can be found 
by the program itself using simple analytic shapes. The grid can consist of multiple blocks, and the program is 
capable of finding its own internal block-to-block boundary surfaces. Volume grid points are found by 
numerically solving the Poisson equations. The Steger & Sorenson (S&S) Right-Hand-Side (RHS) terms (i.e., 
forcing functions) in those equations are of a type which allows the user to choose the desired cell height on a 
read-in boundary, after which the program automatically finds the actual numerical values for the RHS terms 
which yield the desired cell heights. In the process the RHS terms attempt to give local near-orthogonality in the 
region of those same read-in surfaces. The cell heights the user requires may be of any magnitude (limited only 
by the precision of the computer), appropriate for both viscous and inviscid aerodynamic flow modeling. The 
input data is ordinary text, with required formatting. The output grid may be any of three formats, including the 
commonly used PLOT3D 4 formats. 

All the features described above for the original program are preserved in the new program, and a 
significant suite of new features is added. Those new features include: 

• Grid quality is enhanced by re-formulated Steger & Sorenson (S&S) control terms in the Poisson 
Equations. The user may specify arbitrary angles with which lines are to intersect boundaries, rather than 
that specification being limited to 90° everywhere. The treatment of sharp comers which transverse 
boundary surfaces (e.g., a grid wrapping around an airplane fuselage which has a strake) is improved 
using this capability. 

• Another improvement to grid quality is the addition of Thomas & Middlecoff (T&M) clustering terms for 
cases where all six faces of a block are read-in, as found in GRIDGEN and 3DMAGGS. The user can 
choose either the Steger & Sorenson terms (as in the original code and improved as described above), the 
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Thomas & Middlecoff type terms, or a blending between the two which gives good cell-size and 
skewness control at both the boundaries and the interior. 

• Grid quality is evaluated by computing and printing maxima, minima, medians, and averages of cell 
heights and non-orthogonality, at boundaries and in the interiors of the blocks of the finished grid. 

• Initialization is improved by Trans-Finite Interpolation 5 (for cases with six fixed boundary surfaces). In 
some cases grids initialized thusly can serve as the final grid, in others this improved initialization speeds 
convergence. 

• Erlich's Ad Hoc Method for computing locally optimum relaxation parameters is available lor the code's 
SOR solver. This also can accelerate convergence 

• When installed on CRAY computers the code is vectorized in all three coordinate directions, allowing the 
longest possible vector length in each block. This, too, accelerates convergence. 

• The grid generation iteration schedule can be divided into parts. Parameters which effect convergence 
(such as relaxation rates), as well as the type of clustering terms used and their associated decay rates, are 
adjustable with each part. Intermediate solutions and restart files can be written after each part. Thus, in 
practical operation, as much can sometimes be accomplished in one run with this program as in multiple 
runs with other grid generators. 

• An input filter called PREGRAPE/AL, taken from 3DMAGGS, is supplied as a companion program. It 
inputs the output from the GRIDGEN 6 code, which contains blocking strategy and surface grids, and 
turns that into input for 3DGRAPE/AL. 

• Required cell heights and skewness at read-in surfaces can be specified by the user at each point from a 
file. One possible application of the ability to specify the skewness at each point is in hypersonic flow 
where in the exit plane near the shock the flow is aligned with the shock and so should be the grid, while 
in that same face near the body the grid should be aligned with the body surface, and the angles are 
blended in-between. 

• A complete grid generated elsewhere can be read-in, and the elliptic solver can be run a few steps to 
smooth the grid. 

• A Graphical User Interface, coded in FORTRAN-77 and calling the IRIS Graphics Library, allows the 
user to watch selected grid surfaces while the grid solver is iterating. A full suite of transforms and other 
features is included. 

The code is written in FORTRAN-77. It can be installed as an ordinary batch program, and in that form it 
should run on almost any computer. Alternately, on a Silicon Graphics Inc. (SGI) workstation it can be installed 
along with its graphical user interface. The GUI is also written in FORTRAN-77, and calls functions in the IRIS 
Graphics Library. For compiling on a CRAY supercomputer there is a vectorized batch version. 


THEORETICAL DEVELOPMENT OF THE POISSON EQUATIONS IN PHYSICAL SPACE 


The original 3DGRAPE program, the 3DMAGGS program, and the new 3DGRAPE/AL program all 
generate grids by iteratively solving the Poisson Equations in three-dimensions. A mapping is thus found 

between the computational coordinates and the physical coordinates X,Y,Z. The equations are typically 
given in the computational space as 


^x+^yy + ^z = P(^n,0 

(la) 

Uxx 4yy 4zz 0 

(lb) 

Cxx+Cyy+^zz = R(^4,0 

(lc) 
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However, it is natural to apply them in the physical space. It is natural to specify the grid boundary 

conditions by giving X,Y,Z at fixed values of rather than to give values of £,,r|,C at fixed values of X,Y,Z. 
The transformation of Eqs. 1 to physical space proceeds as follows. Clearly, we must have 


% 

h 

C 

To effect this transformation we must also have 

X 

y 

z 


5(x.y.z) 

(2a) 

n(x,y,z) 

(2b) 

C(x,y,z) 

(2c) 

x(£,r|.0 

(3a) 

y(^.r|.C) 

(3b) 

z&ti,© 

(3c) 


Differentiating Eqs. 2 and applying the chain rule gives 




^x Sy ^ 

drj 

= 

Tlx T"|y 

dC 


Cx Cy ^ 


Likewise, differentiating Eqs. 3 and applying the chain rule gives 


dx 

dy 

dz 


dx 


^ X C 

d£ 

dy 

= 

y^ yti yc 

dr| 

dz 


z^ Zr\ ^ 

dC 


(4) 


(5) 


We designate the 3 x 3 matrix in Eq. 5 as M, assume that its inverse exists, and pre-multiply both sides of Eq. 5 
byM' 1 - This gives 



dx 


di; 

M' 1 

dy 

= 

dr) 


dz 




( 6 ) 


Substituting from Eq. 6 into Eq. 4 gives 
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( 7 ) 



dx 


£x £z 

dx 

M' 1 

dy 

= 

rixlly hz 

dy 


dz 


Cx Cy Cz 

dz 


We know that, in general, if 
and if B*^ exists then 
Therefore it must be true that 
Pre-multiplying by B gives 

Applying Eqs. 8 to Eq. 7 gives 


M' 1 = 


Av = Bv 
B 1 A v = v 
B 1 A = I 
A = B 

5 y & 

rixlly Tj z 


(8a) 

(8b) 

(8c) 

(8d) 


( 9 ) 


Cx Cy Cz 

For this to be useful, we must find M~ l . It is known, in general, that 

A -i _ Adj(A) 

Det(A) 

Where Adj(A) is the adjoint of A and Det(A) is the determinant of A. The adjoint of A is a matrix having as each 
element the corresponding cofactor of A. Thus, from Eq. 9, we have 


( 10 ) 


^y Z,z 

r|x ^ly ru 

Cx Cy Cz 


YllYl 2 Yl 3 

Y21Y22Y23/J 

Y3lY32Y33 


( 11 ) 


where Yij is the ij-th cofactor of M and J is the determinant of M. By inspection of Eq. 1 1 we see that 


= Y11/J 

(12a) 

II 

(12b) 
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(12c) 


^z “ Yn/J 

Tlx = Y21/J (12d) 

T| y = Y22/J (12e) 

riz = Y23/J 020 

Cx = Y3l/J ( |2 g) 


Cy = Y32/J 02h) 

Cz = Y33/J (120 

Completion of the derivation of the transformed Poisson equations requires further differentiating the 
metrics in Eqs. 12, substituting them into Eqs. 1, and collecting terms. This process is simple calculus, but very 
lengthy and beyond the scope of this paper. The result is 


where: 


and 


«ll r l£ + CX 2 2^r|r) + OC 33 r^ 

+ 2 ( 0C 12 r^r| + 0Ci 3 r^ + (13a) 

= -J 2 (Pr^ + Qr^i + Rr£) 


X 

r=| y 


“.fJjn.Yn, 


(13b) 


(13c) 


THEORETICAL DEVELOPMENT OF IMPROVED S&S RHS TERMS 

The distribution of points in the grid results primarily from the influence of the Right-Hand Side (RHS) 
terms, or forcing functions. We are free to choose them as we please. In both new and old programs they are: 

= Pidl.Qe-^ + P 2 (Tl,0e-«— 0 

+ P 3 (^,Qe- aT i + P^^Oe- 3 ^™-^) (i4a) 

+ P 5 (^,r|)e-^ + P 6 (^,r|)e- 3 (^ax-0 
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0(5,11,0 = Qi(T|,Qe^ + QsCTl.Oe^S— -0 

+ Q 3 (tOe-«n + Q 4 (^Qe^’w- , i) 
+ Q5(^Tl)e' a? + Q 6 (^ri)e-^™-0 


(14b) 


R(5,T1.0 = Ridl-O®* 5 + Rj(tl,Oe*S--® 

+ ^3(^»O e ' aT1 + R 4 (^Qe- a (w- T i) 

+ R 5 (^,Ti)e-aC + R 6 (^,Ti)e-a(Ua X -Q 


(14c) 


Clearly, these RHS terms P,Q,R are simply superpositions of other terms Pn.Qn^n for l<n<6, multiplied 
by exponentials which are at their maximum value, one, at the boundary surfaces and which decay with distance 
into the interior of the block. The positive constant "a" in Eqs. 14 is set by the user, and determines the rate of 
exponential decay in the size and influence of the RHS terms. 

At this point we must introduce a nomenclature for the face numbers. It is seen in Table I. By examining 
that nomenclature we see that at each of the boundaries the terms in P,Q,R having their subscripts equal to the 
face number are non-zero, and the other terms in P,Q,R approach zero due to the behavior of their exponential 
factors. At face 3, for example, Eqs. 14 reduce to: 


p&ti.O-Pj&O 

05a) 

Q(5.ii,0 = 0*5.0 

05b) 

R&ti,5) = R 3 (5,0 

(15c) 

at face n by considering each face in turn. 

At each point on each face 


we: 

• Assume that the Poisson Equations, Eqs. 13, are satisfied. 

• Find values for all first and second partial derivatives required by Eqs. 13. 

• Eqs. 13 reduce to a 3 x 3 set of linear equations in the three unknowns Pn-On-^n Solve them. 

Having found all the P n ,Qn»Rn, for l^n<6, we can calculate P,Q,R at all points in the grid from Eqs. 14. 

However, finding values for all first and second partial derivatives at each face is not trivial. To further 
illustrate this we must restrict our attention to a particular face. We choose face 3 to illustrate. On face 3 the 

derivatives r^, r£, r^, r^T, and r^ can be found by differencing known boundary face points. The derivatives 

r^ are found by differencing the grid solution at the current time step, as described on page 78 of Ref. 2. If we 

could find derivatives we could then difference them to find derivatives and . 

We find derivatives r^ by adding additional equations which embody the user’s requirements on cell 
height and skewness. In the old 3DGRAPE method we added the three equations 


>-W 

iT Tt 

• 

J 4 

II 

o 

06a) 

? n • ? c = 0 

(16b) 
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(16c) 


^*tj • I*r) — 

As seen in Table 1, E, and C, vary over face 3, and r| varies along lines intersecting the face. Thus Eqs. 16a and 
16b require orthogonality between the lines intersecting the face and the coordinate lines running over the face. 
Eq. 16c requires that the cell height on the surface be the positive constant S. 

It is at this point that the old 3DGRAPE method and the new 3DGRAPE/AL method differ. In the new 
method we realize that when making grids about real-world configurations, with singularities and slope 
discontinuities, it is sometimes necessary to have grid cells which are skewed in a specified way. Lacking this 
ability, an inconsistency can develop which can either cause the elliptic solver to not converge, or result in an 
unsuitable grid. And so Eqs. 16 are replaced by 

. IV, = |r^| . [frj COS 01 ( 17 a) 

?r| . = |rrj • |rc] cos0 2 ( 17 b) 

ft) • r n = S 2 07c) 

where 0j is the angle between the coordinate line intersecting face 3 and the line of varying E, on face 3, and 02 is 

the angle between the coordinate line intersecting face 3 and the line of varying £ on face 3. For 0 \ and 02 equal 

to 90’, Eqs. 17 reduce to Eqs. 16. 

We now proceed to solve Eqs. 17 for r n. Expanding, we have 

Xf;X n + y 5 y n + z f z n = c, 
x ; x ii + y?yn + z c z„ = c 2 

x 3 + yn + A = s 2 

where 

C( = |r^| S cos 0( 
c 2 = |r^| S cos 0 2 

Cj and c 2 are constants because 0 1 , 0 2 , S, and the points on face 3 are user-defined inputs. Equations 18 are 

-* 

three equations in the three unknowns xn^nZn which are the elements of r n. But because Eq. 18c is quadratic, 
solving this set of equations is not straightforward. We will make an assumption about one of the unknowns and 
solve, make that assumption about another of the unknowns and solve, and then make that assumption about the 
last of the unknowns and solve. We will then select the answer which is "best." 

The first assumption we make is that x^ is a constant. Terms involving xq in Eqs. 18a and 18b are 
brought to the right side of the equations, and then the equations are solved, yielding 


y n - x r|Y22 / Yl2 + kl 

(19a) 

z^ = x„Y 32 / Yi 2 + k 2 

(19b) 


where 


(18a) 

(18b) 

(18c) 
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K1 and k2 are constants. Then from Eq. 18c 


kj = 

k 2 = 


ClZ C - c 2 z^ 

I y^ 

c 2 y^ - Ciy c 

-In 


Xn = 


-b ± Vb 2 - 4 ac 
2 a 


(20) 


where 


a = l + (y 22 /y 12) 2 + (Y32/Y12) 2 

b = — (kiYi 2 + k 2 Y 3 2 ) 

112 

c = k 2 + k 2 - S 2 

The second assumption is that y^ is a constant. Terms involving yq in Eqs. 18a and 18b are brought to 
the right side of the equations, and then the equations are solved, yielding 


= yiVYi2 / 722 + ki 
Zt! = ytl732 / 722 + ^2 


(21a) 

(21b) 


where 


kj = 

k 2 = 


CiZ C - c 2 z^ 

Y22 

c 2 x^ - CjX^ 


Then from Eq. 1 8c 


yri = 


Y22 

-b ± Vb 2 - 4 ac 
2 a 


( 22 ) 


where 


a = 1 + (Y12/Y22) 2 + (Y32/Y22) 2 

b = -^-{kiYi 2 + k 2 Y 32 ) 

Y22 


and c is the same as above, in Eq. 20. 
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The third assumption is that zq is a constant. Terms involving zq in Eqs. 18a and 18b are brought to the 
right side of the equations, and then the equations are solved, yielding 


where 


Then from Eq. 17c 


where 


= ZnYu / 732 + ki 

(23a) 

yt| = z x\j22 / J}2 + k 2 

(23b) 

y _ Ciy c - c 2 y^ 

' -Y32 

_ c 2 x^ - C,x c 

2 "732 


-b ± Vb 2 - 4ac 
^ 2a 

(24) 

a = 1 + (Y12/Y32) 2 + (722/732 f 
b - rr - (kiY 12 + ^722) 



and c is the same as above, in Eq. 20. 

In general, none of these three assumptions is strictly correct. However, it usually turns out that at least 
one of them is close enough to correct for this method to generate suitable grids. It was said that we would 
choose whichever of these three solutions was "best." However, Eqs. 20, 22, and 24 each include an ambiguous 
sign from a square-root operation. Therefore, we actually have six solutions to choose from. Using each of the 
six solutions we compute the Jacobian. If the coordinates in the block are right-handed (with the "handedness" 
being a user-defined input) we choose the solution which yields the largest positive Jacobian. If the coordinates 
in the block are left-handed we choose the solution which yields the largest negative Jacobian. The logic behind 
choosing based upon the Jacobian is that Jacobians, as defined above, having large absolute values seem to be 

present in grids which are more orthogonal, and, conversely, Jacobians having small absolute values seem to be 

-► 

present in grids which are highly skewed. Thus the elements of fit are found. 

The foregoing is the analysis for face 3. The analysis for face 4 appears identical, differing only in some 
of the difference formulas. The analyses for faces 1 , 2, 5, and 6 follow in a straightforward manner from the 
foregoing example. 

This formulation for the S&S RHS terms requires a lot of computation but most of it is done only once, at 
the start of the iteration schedule. It was said, above, that having all values for the derivatives at the face, those 
derivatives are substituted into Eqs. 1 3, yielding a 3 x 3 set of linear equations in the three unknowns P n ,Q n J*n 

Their solution shows Pn.Qn*Pn to be linear functions of the second derivatives r^ which are found by 
differencing at each time step. The coefficients in those linear functions are fixed for all computational time. 
Therefore, the only computation necessary to find the RHS terms in each iteration is to re-evaluate r^ , re- 
evaluate the linear functions to get PrnQn>Pn at each face, and then use Eqs. 14 to re-compute the P,Q,R at every 
point in the grid. 
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The effectiveness of this method is seen in Figure 1. When wrapping a grid around a sharp edge it is 
necessary to cause the lines intersecting the surface near the edge to bend toward the edge for best results. The 
ultimate example of wrapping a grid around a sharp edge is to wrap it around the edge of a flat plate. Figure 1 
shows a wing with a zero-thickness extension in the spanwise direction, and a C-H type grid around it. Thus, it is 
necessary to wrap a C-type grid around the leading-edge of that wing and its flat-plate extension. This would not 
have been possible with the old type RHS terms. 


THOMAS AND MIDDLECOFF CLUSTERING TERMS 

When making grids in regions where all six faces of the computational cube are fixed it is sometimes 
advantageous to use clustering functions where the spacing normal to a face is determined by the spacing on the 
side walls. The Thomas and Middlecoff 7 clustering terms are included here for that purpose. However, the 
Thomas and Middlecoff clustering terms 


where 


P = <D(V^ . V£) 

(25a) 

Q = 'P(Vt) . Vr|) 

(25b) 

R = Q(VC . VO 

(25c) 

• % 

(23d) 

< 

n 

j* J 4 
• 

j* i" 

(25e) 

a-*-!? 



r c . r c 


are given in the computational space, and to be useful here they must be converted to physical space. Applying 
the definition of the V operator, illustrated by 


Vq = y + E, y k + U a® 

A /N A 

where j, k, and 1 are the unit normal vectors, and reducing, gives 

P = + & (27a) 

Q = vKhx + n? + ni) ( 27 b) 

R = £2(£* + Cy + C) (270 


Substituting the metrics shown in Eqs. 12 into Eqs. 27, and expanding and re-grouping, gives 
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(28a) 


P = <&[(?„ . ?„) (r c . r c ) - (r„ . \p 2 
Q = V[(^ • ? 5 ) (r c . r ; ) - (r 5 . ?;)]/J 2 ww 

R = £2[(?|= . r 5 ) (r„ . r n ) - . r n )]/J 2 am 

These RHS terms generate good grids in many applications. An exception is the situation where the 
opposing side boundaries, from which the T&M terms are calculated, have very different clustering 
characteristics. In these cases instabilities in the Poisson solver can result. 

It was found in the development of 3DMAGGS that S&S clustering terms tend to give the most-nearly- 
orthogonal grids near boundaries, while T&M clustering terms give the best clustering in the interior of the 
blocks. And so a blending between the two kinds of RHS terms was developed, and is included in 
3DGRAPE/AL. 


OPTIMUM RELAXATION PARAMETER 


3DGRAPE/AL solves the 3-D Poisson equations using Point Successive Over Relaxation (PSOR). In 
PSOR there is a relaxation parameter, £2, which determines the rate of convergence and stability of the method. 

In the old program the £2 was fixed for all computational time. That option is still available in the new code as 
well. However, the new code also has an algorithm to compute an optimum relaxation parameter at every point in 

the grid using the method of Erlich.^ 

That method requires the equations being solved, here Eq. 13a, to be represented as a difference equation 
of the following form: 


ao?j,k,i + aiFj+uj + ^*+1,1 + ^^kj+i 

+ a 4?j-l,k,l + a 5?j,k-l,l + a 6rj,k,l-l = bj,k,l 


(29) 


Applying standard central differences to all first and second partial derivatives in Eq. 13a, and collecting terms, 
we arrive at the form of Eq. 29, where 


a. = .2l3i,+ 3k + 0333 


m w w 


a l - 


a 2 - 


a 3 = 


an , 

J 2 P 


2 A£ 

a 22 , 

J 2 Q 

(An ) 2 

2 Ar| 

a 33 , 

J 2 R 

(AC ) 2 

2 A £ 


(30a) 


(30b) 


(30c) 


(30d) 
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(30e) 


„ _ an J 2 P 

4 (A^) 2 2A^ 

„ _ CX 22 J2 Q 

(Ar| f 2Ati 
a - a 33 J 2 R 

6 " (AC) 2 ' 2Ai; 

The complex eigenvalues of Eq. 29 at each point, ignoring wave numbers above 1, are 


(300 


(30g) 


u = (X r + Hi = /aTaJcos^^-- + fai ajcos n 

^ Jmax ' * K max T 1 (31 ) 

+ Vascos ^ — ) 

*max * 1 

where p r and pi are the real and imaginary parts of p, respectively. It is required that |pt| ( 1- 

Continuing with Erlich's method, as formulated by Steinbrenner, Chawner, and Fouts on pages 6-6 and 6- 


7 of Ref. 6 (with typographical errors corrected), we let 

A = |1 2 + |1 2 (32a) 

B = (32b) 

C = A 2 - B 2 (32c) 

D = A 2 - B (32d) 

E = VC + D 2 (32e) 

F = VC (320 


Then 

- (3D + E) F V&T) - (3D - E) F Ve+D + A 2 + 3B 2 - 4 A 2 B 

®~ AT) <3 

and the relaxation parameter to is 

f -(co - Vo? + 4to)/2 if D>0\ 

W |-(oi + Vo? + 4a>)/2 ifD(0( 

This method can reduce the number of iterations required to achieve convergence. The to so computed 
can sometimes be a little too large, and so cause instabilities. Therefore, in the code, they are multiplied by a 
limiting factor. The default value of this factor is 0.7, but a value of 0.6 was found to be necessary in one of the 
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sample cases. These w are dependent on the grid at its current time step, and so they are re-calculated each time 
step. As can be inferred from the above, computing them requires a significant amount of computer time, so the 
code has an option wherein they are re-calculated every n time steps. 


CURRENT AND FUTURE WORK 


Extensions and improvements to 3DGRAPE/AL planned for the near future, or currently in progress, 

include: 


• Adding cylindrical topology -- solving in p,9 ,x — to improve solutions about cylindrical axes, just as use 
of spherical topology is offered to improve solutions about polar axes. 

• Write a translator to translate input for the old program 3DGRAPE into input for the new program 
3DGRAPE/AL. This will make datasets used in the old program workable in the new program. Also, the 
latest release of GRIDGEN, Version 9, can generate input to the earlier 3DGRAPE program, and so this 
translator will provide an alternate pathway between the two codes. 

• We will try again to build multigrid into the Poisson solver. We tried once before, but did not succeed. 
We know that in grids produced by this code the cell size increases in an approximately exponential 
fashion with distance from a boundary, but we cannot predict exactly what that rate of increase in cell 
height will be. Therefore, although we know the desired distance to the 1st node away from the 
boundary, we cannot say, exactly, what should be the distances to the 2nd, 4th, 8th, 16th, etc., nodes from 
the boundary. But we need those distances to operate at the various multigrid levels. And so we 
estimate. But those estimates are not accurate. Thus, there are inconsistencies between the equations 
being solved at the different multigrid levels, and therefore it does not converge. But we are determined 
to build in a multigrid capability if at all possible. 

• Optimize the code to run faster on SGI machines, such as ONYX. Portability to other platforms will be 
preserved. 

• Write a GUI for PREGRAPE/AL. 

• Put in a new boundary condition wherein the angles of lines intersecting a fixed internal boundary is 
mimicked across that boundary. 

• Make a version of the code which uses PVM to run the code in parallel on multiple platforms. 
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Table I. Face numbers. 


Face no. n : 
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K 

n 
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1 
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varying 

varying 
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varying 
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jmax 
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Cmax 

varying 

varying 
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varying 

i 

vaiying 

varying 
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varying 

4 
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varying 

varying 

Ornax 

varying 
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varying 

varying 

i 

varying 

vaiying 

0. 

6 

varying 

varying 

Imax 

varying 

varying 

Cm ax 
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Figure la. Wing With Flat-Plate Extension in Both 
Span wise and Downstream Directions. 



Figure lb. View From Outboard End Showing Portions of Two 
Interior Grid Surfaces Wrapping Around Wing and Extension 
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Figure lc. Close-Up of Grid Surface Wrapping Around 
Leading Edge of Flat-Plate Extension 
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